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Abstract 



A new pseudoperturbative ( artificial in nature) methodical pro- 
posal [15] is used to solve for Schrodinger equation with a class of 
phenomenologically useful and methodically challenging anharmonic 
oscillator potentials V(q) = a Q q 2 + aq A . The effect of the [4,5] Pade 
approximant on the leading eigenenergy term is studied. Comparison 
with results from numerical ( exact) and several eligible ( approxima- 
tion) methods is made. 
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1 Introduction 



Quartic anharmonic interactions continue to remain a focus of attention. 
Their Hamiltonian 

H=£- + a r 2 + ar A (1) 

forms one of the most popular theoretical laboratories for examining the va- 
lidity of various approximation techniques and represents a nontrivial physics. 
Interest in this model Hamiltonian arises in quantum field theory and molec- 
ular physics [1-6]. 

Although enormous progress has been made over the years in our under- 
standing of this Hamiltonian, questions of delicate nature inevitably arise in 
the process. The hardest amongst often relate to the existence of the assumed 
small expansion parameter and the universality of an adequately attendant 
powerful approximation. The implementation of Rayleigh- Schrodinger per- 
turbation theory, or even naive perturbation series, expresses the eigenvalues 
as a formal power series in a which is quite often divergent, or at best asymp- 
totic, for every a^O. One has therefore to sum up such series [7-10]. Hence, 
apparently artificial perturbation recipes have been devised and shown to be 
ways to make progress [2,3,11-16]. Without being exhaustive, several eligible 
methods have been used to calculate the eigenvalues and eigenfunctions for 
Hamiltonian (1). Long lists of these could be found in Ref.s[2, 3, 8-10, 13,17- 
19]. 

In this paper we introduce, in section 2, a new analytical ( or, preferably, 
semianalytical) perturbation method for solving Schrodinger equation. The 
construction of which starts with the time-independent one-dimensional form 
of Schrodinger equation, in h — m — 1 units, 
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2dq 2 



1(1 + 1) 
2q 2 



+ V(q) V nr ,l(q) = E nrtl V nril (q), 



(2) 



where / is some quantum number and n r counts the nodal zeros in ^n r ,i(q)- 
The symmetry of an attendant problem obviously manifests the admissi- 
bility of the quantum number I: In one-dimension (ID), / specifies parity, 
(— with the permissible values -1 and/or ( even and/or odd parity, 
respectively) where q = x G (—00,00). For two-dimensional (2D) cylindri- 
cally symmetric Schrodinger equation one sets I — \m\ — 1/2, where m is the 
magnetic quantum number and q = (x 2 +y 2 ) 1 ^ 2 £ (0, 00). Finally, for three- 
dimensional (3D) spherically symmetric Schrodinger equation, / denotes the 
angular momentum quantum number with q = (x 2 + y 2 + z 2 ) 1 ^ 2 E (0, 00). 

We shall focus our attention, in section 3, on ID and 3D problems and 
consider, for the sake of diversity; (i) 3D anharmonic oscillators V(r) = 
r 2 /2 + r 4 /2 with n r = and I = 0,1,2,5,10,50, (ii) 3D ground state, or 
equivalently ID first excited ( odd-parity) state, for anharmonic oscillators 
V(q) = <? 2 /2 + a? 4 over a wide range of anharmonicities ( i.e.; a = 0.002 
to a = 20000), and (iii) 3D single- well anharmonic oscillator ground state, 
or equivalently ID double-well anharmonic oscillator first excited state, for 
V(q) = -aq 2 /2+q A /2 at various well depths ( i.e.; a = 1, 5, 10, 15, 25, 50, 100). 
For the sake of comparison, we use results from exact numerical methods 
reported in [2,5], the best estimation of the phase- integral method (PIM) [5], 
an open perturbation technique [2], and a perturbative- variational method 
(PVM) [6] . Section 4 is reserved for concluding remarks. 
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2 The Method 



Our methodical proposal uses 1/1 as a perturbation expansion parameter, 
where I = I — f3 and (5 is a suitable shift mainly introduced to avoid the trivial 
case I = 0. Hence, hereafter, it will be referred to as the pseudoperturbative 
( artificial in nature) shifted-/ expansion technique (PSLET). Equation (2) 
thus becomes 



IA 2 

2dq 



2 +V(q)\V nr:l (q) = E nr ^ nrtl ( q ), (3) 



Vfa) = P + ( ^ + g +W + 1) + gv (g ). (4) 



Herein, it should be noted that Q is a constant that scales the potential V(q) 
at large - / limit and is set, for any specific choice of / and n r , equal to I 2 at 
the end of the calculations [11,16]. And, (3 is to be determined in the sequel. 
PSLET procedure begins with shifting the origin of the coordinate through 



x = l 1/2 (q-qo)/q , (5) 



where q Q is currently an arbitrary point to perform Taylor expansions about, 
with its particular value to be determined. Expansions about this point, 
x = (i.e. q = q ), yield 

T To 
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vw.))-f(^)^. 

Obviously, the expansions in (6) and (7) center the problem at an arbitrary 
point q Q and the derivatives, in effect, contain information not only at q a but 
also at any point on g-axis, in accordance with Taylor's theorem. Also it 
should be mentioned here that the scaled coordinate, equation (5), has no 
effect on the energy eigenvalues, which are coordinate - independent. It just 
facilitates the calculations of both the energy eigenvalues and eigenfunctions. 
It is also convenient to expand E as 

oo 

E n r ,l = E E ( S~ n . (8) 
n=-2 

Equation (3) thus becomes 

V^x) = S-E nr ^ nr:l (x) : (9) 



Id 2 a 2 ~ 

-2& + T v <*<» 



with 



I 



V(x(q)) = q 2 J 



1 V(qo) 
2q 2 Q 



+ V-' 2 



-x + 



Q 



+ 



V"{<lo)<l oX 2 

2Q 



-x" + 



+ (2/3 + !)£(-!)' 



(n + 1) 



-x 



+ ? 2 E 



n=3 



n=l 

'd n V{q )\ (q xY 



dq™ ) n\Q 



j-(n-2)/2 



+ {3({3+l)±(-l) ni -^xH-^ 2 +W±^, 



(10) 



n=0 
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where the prime of V(q Q ) denotes derivative with respect to q Q . Equation (9) 
is exactly of the type of Schrodinger equation for one - dimensional anhar- 
monic oscillator 



x n j x ) = vi.,14 (ii) 



where P(x) is a perturbation - like term and e Q is a constant. A simple 
comparison between Eqs.(9), (10) and (11) implies 



e Q = I 



I qlv{q ) 

2 



Q 



2(3+1 P(J3 + 1) 
2 + 21 



(12) 



^n r — ^ 



+ 



I, qjnsA 

2 Q 
g±l) + A (0 



2/3+1 , 1 ' 
2 + K + 

oo 

+ E AS- 1 )/-, 

n=2 



(13) 



and 



A - o 2 V ft (n) H n+1 ) 



n=-2 



(14) 



Equations (13) and (14) yield 



E { ~ 2) - — 



1 + V(? ) 



2g 2 Q 



(15) 



i^ -1 ) = 1 

Ho 



— 2 — + K + 



(16) 
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F (o) = i r w+i) , wo 

n r ,l 2 9 «r 

yo L z 



(17) 



ra > 1. 



(18) 



Here g G is chosen to minimize E ri 




dE, 



;(-2) 
'n r ,l 



= and 




> 0. 



(19) 



dg 



Hereby, V(q) is assumed to be well behaved so that ^ has a minimum q Q 



and there are well - defined bound - states. Equation (19) in turn gives, with 



Consequently, the second term in Eq.(10) vanishes and the first term adds 
a constant to the energy eigenvalues. It should be noted that energy term 
l 2 E^~ r f has its counterpart in classical mechanics. It corresponds roughly to 
the energy of a classical particle with angular momentum L z —l executing cir- 
cular motion of radius q Q in the potential V(q a ). This term thus identifies the 
leading - order approximation, to all eigenvalues, as a classical approxima- 
tion and the higher - order corrections as quantum fluctuations around the 
minimum q , organized in inverse powers of I. The next leading correction to 
the energy series, lE^^, consists of a constant term and the exact eigenval- 
ues of the unperturbed harmonic oscillator potential w 2 x 2 /2. The shifting 




(20) 
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parameter f3 is determined by choosing lE^^—O. This choice is physically 
motivated. It requires not only the agreements between PSLET eigenvalues 
and the exact known ones for the harmonic oscillator and Coulomb potentials 
but also between the eigenfunctions. Hence 

(21) 



- + K + -)w 



where 



w 



qoV"(qo) 
V(Qo) 



3 + 



(22) 



Then equation (10) reduces to 



f V(x(q)) = qll 



1 V(Qq 

H Q 



(23) 



n=0 



where 



/J^^^Vi^, (24) 



v^{x) = -(2/3 + l)x-2x 3 +^ V {qo) - 3 



6Q 



x 



(25) 



and for n > 2 



; _l)«( 2/ 3 + l)^±^-x n + (-i)»^±ll ( n - l)^"- 2 ) 
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+ (-1) 



n 



(71 + 3) 

2 



+ 



g (n+4) d n+2 V(q ) 
Q{n + 2)\ dqZ+ 2 



x 



,n+2 



(26) 



Equation (9) thus becomes 




+ £yn)j-»/2 * np>| ( x ) = 



(27) 



Up to this point, one would conclude that the above procedure is nothing 
but an imitation of the eminent shifted large-N expansion (SLNT) [12,14,16,20- 
22]. However, because of the limited capability of SLNT in handling large- 
order corrections via the standard Rayleigh-Schrodinger perturbation theory, 
only low-order corrections have been reported, sacrificing in effect its precise- 
ness. Therefore, one should seek for an alternative and proceed by setting 
the nodeless, n r = 0, wave functions as 



In turn, equation (27) readily transforms into the following Riccati equation 
[2,3, and references therein]: 



*o,i(s(5)) = exp(U 0>l (x)). 



(28) 



-\U"{x) + U'{x)U'{x)\ + ]T v {n \x)l- n/2 



n=0 



1 f 13(13 + 1) 




) 



CO 

+ £ \t 1] l 



—n 



(29) 



n=2 
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Hereafter, we shall use U(x) instead of Uo,i(x) for simplicity, and the prime 
of U (x) denotes derivative with respect to x. It is evident that this equation 
admits solution of the form 



where 



u\x) = jr u {n \x)r n ' 2 + jr G (n) (^)r (n+1)/2 , 



n=0 



n=0 



n+1 



2m- 1 



m=0 



.X 



n+1 



; A),n — o, 



G (n) (l) = £ C m,nX 



2m 



m=0 



Substituting equations (30) - (32) into equation (29) implies 



(30) 



(31) 



(32) 



1 oo 

n=0 

oo oo 



jj(n) i~n/2 + ^(n) j-( n +l)/2 
1 oo oo 



n=0p=0 



+ g w (n) r -»/2 = l| r ^^ + ^oA + g^n-l) r - rlj (33) 
n=0 t \ 2 / n=2 

where primes of U^ n \x) and G^rr) denote derivatives with respect to x. 
Equating the coefficients of the same powers of I and x, respectively, ( of 
course the other way around would work equally well) one obtains 



I[/(0)' _i u (o) u (o) +v (o) =0 
2 2 



(34) 
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U<®\x) = D lt0 ; D lfi = 
and integration over x yields 

U {0 \x) = -wx. 



-w 



Similarly, 



![[/(!)' + G (0)'] _ ^(0)^(1) _ ^(0)^(0) + 

2 



U<V(x) = 0, 



G^(x) = C 0fi + C lfi x\ 



C - Bl 



C , o,o = -(Ci, + 2/3 + 1), 
w 
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_I[[/( 2 )' + G {1) '] - - X] C/ (n) C/ (2_n) - -G (0) G (0) 

2 2 n=0 2 

- £ [/W^ 1 "") + w ( 2 ) = ^-til + A (o) > (43) 

n=0 2 



f/( 2 )(:r) = D li2 i + V- ( 44 ) 



G (1 \x) = 0, (45) 



2 



^2,2 = -(-^ - Sa) (46) 



A,2 = ^(^2,2 + C ,oC li0 - ?(2/J + 1)), (47) 



5 g D 6 rfV(g ) 

52 " 2 + 24g^T' (48) 



Ai 0) = -^(A,2 + C 2 ). (49) 



and so on. Thus, one can calculate the energy eigenvalue and the eigenfunc- 
tions from the knowledge of C m ,n and D mn in a hierarchical manner. Nev- 
ertheless, the procedure just described is suitable for systematic calculations 
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using software packages (such as MATHEMATICA, MAPLE, or REDUCE) 
to determine the energy eigenvalue and eigenfunction corrections up to any 
order of the pseudoperturbation series. 

Although the energy series, Eq.(8), could appear divergent, or, at best, 
asymptotic for small /, one can still calculate the eigenenergies to a very good 
accuracy by forming the sophisticated [N,M+1] Pade' approximation 

P^ + \l/T) = (P + Pi/l+ ■■■ + Pm/1 M )/(1 + qi/l + ■■■ + q N /l N ) 

to the energy series [23]. The energy series, Eq.(8), is calculated up to E^J /I s 
by 

Eo,i = PE^ + E$ + ■■■ + E$/P + 0(1/P), (50) 

and with the P|(l/Z) Pade' approximant it becomes 

E , l [A,5]=PE i - 2) + P!(l/l). (51) 



3 Quartic anharmonic interactions 

Let us consider the phenomenologically useful and methodically challenging 
quartic anharmonic interactions 

V(q) = a q 2 + aq 4 (52) 
of Hamiltonian (1). Equation (22) then reads 
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8a q + 2Aatql 
2a q + Aaql ' 



(53) 



and Eq.(20) yields 



8a q + 2Aaql 
\ 2a q o + Aaql 



= q 2 ^2a + Aaq 2 



(54) 



In the absence of a closed form solution for q Q in (54), one should appeal 
to some software packages ( MAPLE is used here) to resolve this issue. Of 
course there is always more than one root for (54). However, the symmetry 
of the problem in hand along with Eq.(19) would single out one eligible root 
q Q as a minimum of E^' 2 \ Once q Q is determined the coefficients C m ^ n and 
D m , n are obtained in a sequential manner. Consequently, the eigenvalues, 
Eq.(50), and eigenfunctions, Eqs.(30)-(32), are calculated in the same batch 
for each value of a Q , a, and I. 

Our results ( tables 1-3) are obtained from the first eleven terms of our 
energy series (50). Also, the effect of the [4,5] Pade approximant on the 
leading term PE^~ 2 ^ is reported as E[4,5]. In table 1 we list our results along 
with the exact numerical ones and the (best estimated) eigenvalues obtained 
from the fifth-order phase-integral method (PIM) reported by Lakshmanen et 
al. [5]. Obviously, our results compare excellently with the exact numerical 
ones and surpass those from PIM. Whilst the [4,5] Pade approximant had 
no dramatic effect on the energy eigenvalues for I — 0, it had no effect on 
the energy eigenvalues for I > 1. A common feature between PSLET and 
PIM is well pronounced here; the precession of both methods increases as I 
increases. 

Again we proceed with the theoretical laboratory (52) and examine the 
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validity of PSLET over a wide range of anharmonicities for V(q) = q 2 /2 + 
aq 4 . In table 2 we list our results for the three-dimensional (3D) ground 
states energies, or equivalently for the one- dimensional (ID) first excited state 
energies. The results of Bessis and Bessis [2], via an open perturbation recipe, 
and the exact ones [24], using Bargman representation, are also displayed. 
Clearly and satisfactorily, the trend of the exact values of the energies is 
reproduced. 

Finally, we consider the ground state energies of the 3D single-well, or 
equivalently the first excited state energies of the ID double-well, potentials 
V(q) = — ag 2 /2 + g 4 /2. We compare our results ( table 3) with those obtained 
by Saavedra and Buendia [6] via a perturbative- variational method (PVM). 
They are in excellent agreement not only with the PVM but also with the 
hypervirial perturbation method [25], especially for deep wells. 

4 Concluding remarks 

The method (PSLET) just described is conceptually sound. It avoids trou- 
blesome questions such as those pertaining to the nature of small-parameter 
expansions, the trend of convergence to the exact numerical values, the utility 
in calculating the eigenvalues and eigenfunctions (in one batch) to sufficiently 
heigher-orders, and the applicability to a wide rang of potentials. Provided 
that the latter is analytic and give rise to one minimum of E^~ 2 ^ and an 
infinite number of bound states. 

On the computational and practical methodology sides, PSLET comes in 
quite handy and very accurate numerical results are obtained. Nevertheless, 
if greater accuracy is in demand, another suitable criterion for choosing the 
value of the shift f3, reported in [13,25], is also feasible. However, one would 
always be interested, for practical exploratory purposes, in the conventional 
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wisdom of perturbation prescriptions that only a few terms of a "most useful" 
perturbation series reveal the important features of the solution before a state 
of exhaustion is reached. Our method indeed belongs to this category where 
the results of the illustrative challenging examples used bear this out. 

On the other hand, asymptotic wavefunctions emerge in our procedure 
from the knowledge of C m n and D m n to study, for example, electronic tran- 
sitions and multiphoton emission occurring in atomic systems. Such studies 
already lie beyond the scope of our present methodical proposal. 
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Table 1: Eigenvalues from the fifth-order phase-integral method E PIM [5], 
the pseudoperturbative shifted-/ expansion technique Epslet, the effect of 
the [4,5] Pade approximant on our leading energy term E[4,5], and from 
the exact numerical calculations [5] for the three-dimensional anharmonic 
oscillator V(r) = \r 2 + ^r 4 , with n r = 0. 





1=0 


1=1 


1=2 


EpiM 


2.324 83 


4.190 26 


6.242 80 


EpsLET 


2.324 40 


4.190 17 


6.242 78 


E[4,5] 


2.324 41 


4.190 17 


6.242 78 


E eX act 


2.324 41 


4.190 17 


6.242 78 




1=5 


1=10 


/=50 


EpiM 


13.264 459 9 


27.092 492 362 


187.529 708 014 021 


EpsLET 


13.264 458 8 


27.092 492 304 


187.529 708 014 0025 


E[4,5] 


13.264 458 8 


27.092 492 304 


187.529 708 014 0025 


E eX act 


13.264 458 8 


27.092 492 305 


187.529 708 014 003 
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Table 2: Three-dimensional ground state energies or equivalently one- 
dimensional first excited state energies for V(q) = y + aq 4 . Ebb denotes 
Bessis and Bessis results [2] and the exact ones E exact , reported therein, for 
different anharmonicities. 



a 


EpsLET 


£[4,5] 


Ebb 


E eX act 


0.002 


1.507 41940 


1.507 41940 


1.507 4194 


1.507 41939 


0.006 


1.521 80570 


1.521 80570 


1.521 8057 


1.521 80565 


0.01 


1.535 64844 


1.535 64846 


1.535 6483 


1.535 64828 


0.05 


1.653 439 


1.653 439 


1.653 441 


1.653 43601 


0.1 


1.769 512 


1.769 625 


1.769 529 


1.769 50264 


0.3 


2.094 678 


2.094 640 


2.094 795 


2.094 64199 


0.5 


2.324 401 


2.324 407 


2.324 661 


2.324 40635 


0.7 


2.509 16 


2.509 23 


2.509 56 


2.509 22810 


1 


2.737 73 


2.737 91 


2.738 32 


2.737 89227 


2 


3.292 48 


3.292 94 


3.293 50 


3.292 86782 


50 


8.913 21 


8.916 61 


8.917 41 


8.915 09636 


200 


14.056 17 


14.062 53 


14.062 96 


14.059 2268 


1000 


23.966 93 


23.978 93 


23.978 63 


23.972 2061 


8000 


47.880 19 


47.890 95 


47.903 66 


47.890 7687 


20000 


64.972 32 


65.006 64 


65.004 18 


64.986 6757 
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Table 3: Three-dimensional ground state energies or equivalently one- 
dimensional first excited state energies for V(q) = —aq 2 /2 + g 4 /2. E PVM 
represents the results from perturbative- variational method [6]. 



a 


EpvM 




E[4,5] 


1 


2.834 5 


2.835 3 


2.834 4 


5 


-3.250 68 


-3.250 85 


-3.250 84 


10 


-20.633 55 


-25.633 69 


-20.633 50 


15 


-50.841 387 


-50.841 42 


-50.841 42 


25 


-149.219 456 


-149.219 454 


-149.219 454 


50 


-615.020 090 9 


-615.020 091 


-615.020 091 


100 


-2845.867 880 34 


-2485.867 880 337 


-2485.867 880 337 
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